Packing structure of a two-dimensional granular system through the jamming 
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We have performed a novel experiment on granular packs composed of automatically swelling 
particles. By analyzing the Voronoi structure of packs going through the jamming transition, we 
show that the local configuration of a jamming pack is strikingly similar to that of a glass-forming 
liquid, both in terms of their universal area distribution and the process of defect annealing. Fur- 
thermore, we demonstrate that an unambiguous structural signature of the jamming transition can 
be obtained from the pair correlation functions of a pack. Our study provides insights into the 
structural properties of general jamming systems. 
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Nearly 300 years ago, Stephen Hales, a highly gifted 
English scientist and clergyman, designed a clever ex- 
periment to investigate the packing structure of spheri- 
cal particles He put spherical peas into a closed jar 
filled with water. Absorbing the water slowly, the peas 
confined in the jar expanded and deformed under the 
enormous stress of swelling. From the shape of deformed 
peas, Hales deduced the structure of the initial pea pack. 
Though this was probably not Hales' intention, his clas- 
sic experiment provides an excellent way to probe the 
jamming transition of a granular pack. In this paper, we 
report a two dimensional (2D) version of Stephen Hales' 
experiment. From the position of swelling particles in a 
cell of fixed area, we investigated the structure of a gran- 
ular pack as it evolves continuously through the jamming 
transition. 

The jamming transition and the jamming phase dia- 
gram are intriguing concepts for unifying different tran- 
sitions in various systems ranging from granular mat- 
ter to molecular glasses Many interesting features 
have been found in athermal granular systems near the 
jamming transition However, most of these re- 
sults were obtained only in theories and simulations for 
ideal systems of frictionless spherical particles. It is not 
straightforward how such results should be extended to 
experimental granular systems with frictional contacts. 
Majmudar et al. found that the increasing of both the av- 
erage coordination number of particles and the pressure 
of the system agree with the mean-field theory of friction- 
less particles as the jamming transition is approached 
However, other measurements produce different results 
in experiments and simulations I3~12|. For example, 
experiments showed that the ratio of the shear modu- 
lus to the bulk modulus of a pack, G/B, stays constant 
as the pressure of the system approaches zero 13, 11|. 
which contradicts the results of simulations and theories 
with frictionless particles where G/B diminishes with the 
pressure [H, Moreover, very few experiments have di- 



rectly investigated the structural signature of the jam- 
ming transition. Corwin et al. indirectly deduced the 
structural signature of the transition from the statistics 
of the contact forces between particles [TH • Other than 
granular systems, Zhang et al. investigated the structural 
signature of a colloidal system at finite temperature 14 1 . 

Here, we experimentally study the structure of a gran- 
ular pack going through the jamming transition. Instead 
of peas, we used tapioca pearls as our granular material. 
Tapioca pearls are spherical particles made of starch that 
swell up by absorbing water. The average diameter of 
dry pearls is do = 3.3 mm with the polydispersity of 8%, 
which prevents the system from forming large crystalline 
areas. A fully swelled particle can be as large as 1.7 times 
of its original size, which leads to almost 3 times increases 
in packing fraction ((f)) of a 2D system. More impor- 
tantly, the swelling of pearls is slow and uniform, and 
the pearls maintain their spherical shape during swelling. 
The polydispersity of swelled particles becomes slightly 
smaller than that of dry particles. 

For each experiment, a single layer of pearls is laid 
down randomly into a square cell with the side length 
of L = 54.6 cm. The vertical height of the cell is kept 
at 5.454 mm using a spacer and washers, which prevents 
fully swelled pearls from buckling into two layers (Fig.QJi 
inset). Since the jamming transition occurs before pearls 
fully swell, they are not stuck between the two confining 
plates near the jamming point. The cell can hold over 
15,000 particles at the initial packing fraction 4>q > 0.63. 
A subset of the sample (about 10,000 particles) is stud- 
ied in the central area of the cell to avoid any boundary 
effects. To allow the pearls to swell, the whole system in- 
cluding the cell and particles is submerged under water. 
Static friction between particles and the bottom plate is 
also reduced due to the lubrication of water. A high res- 
olution camera is used to take an image of the sample 
every 20 s. Due to the slow swelling of pearls, the system 
is quasi-stationary at this time scale. The acquired im- 
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age s are then analyzed with a particle tracking algorithm 
[15j. The center of particles can be determined with an 
error of 0.1 mm (about 3% of a particle diameter). A 2D 
projection of images is used to obtain </>. More details of 
the experiment can be found in 
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Wc initially prepared the sample in a dense but un- 
jammed state. As the size of particles increases, the sys- 
tem goes through the jamming transition and is eventu- 
ally trapped inside the jammed phase. We constructed a 
Voronoi tessellation from the centers of particles for each 
image with increasing <fi. The Voronoi cell of a particle is 
the polygon consisting of all points closer to the center of 
the particle than to centers of any other particles, which 
reveals local configurations of packings [1| . 

We first studied the area distribution of the Voronoi 
cells, P(s), at different <fi as the system goes through the 
jamming transition. We use a scaling that has been pre- 
viously used in analyzing the structure of glass-forming 
liquids [13] ■ The area of each cell s is shifted by the 
average area of all cells (s) and then is normalized by 
the standard deviation a s , i.e. s — (s — (s))/a s . The 
probability distribution function is then P(s) — <r s ■ P(s) 
(Fig. [1]). The shape of normalized distributions at small 
<t> depends on the initial random configuration of parti- 
cles (Fig. QJ,), which varies for each run. However, as 
(j) increases toward the jamming point, we found that 
all P(s) collapse into a single master curve. For clar- 
ity, the collapsed P(s) is plotted separately in Fig. [Tp 
and also in a log-linear plot in the inset of Fig. [TJd. It 
is surprising that all the packs with (f> ranging contin- 
uously from 0.83 up to 0.90 (which is the upper limit 
of our experiment due to the finite swelling potential of 
particles under compression) follow the same universal 
distribution. From the structural signature of the jam- 
ming transition which we shall show later, we know that 
the system jams at <jf ~ 0.84. Hence, the universal dis- 
tribution emerges only close to and after the jamming 
transition. It should be emphasized that the universal 
distribution for this large range of <fi is non-trivial: the 
jammed phase of a pack is far from stationary. Though 
particles are confined by their neighbors, the shape of 
Voronoi cells is highly mobile as illustrated below. Some 
localized collective rearrangement of particles also exists 
in the jammed phase (16J. 

Quantitatively, wc fitted the universal distribution 
with a shifted gamma distribution (Fig. Hb), 



p(s) = (s- s y 



_ 1 exp(-(s - s )/ 6) 
T(k)8 k 



(1) 



where 9 and k are standard parameters of gamma distri- 
bution and so is a fitting parameter due to the shift in 
the s axis. Th e same distribution has been suggested be- 
fore in 
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1 91 ] . We also compare the gamma distribution 
with a Gaussian distribution (Fig.QJ)), and find that due 
to the asymmetric shape of the distribution, the gamma 
distribution gives a better description of our data. 
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FIG. 1: Distribution of Voronoi cell areas, shifted by the aver- 
age area and scaled by its standard deviation, (a) Distribution 
for a pack with <f> increasing from 0.70 to 0.90. Upper right in- 
set shows a schematic of the experimental setup. The spacer 
and washers keep the particles in one layer. The washers al- 
low water to flow in and out of the otherwise closed cell, (b) 
The collapse of the distributions near and after the jamming 
transition. The universal distribution of Voronoi structure 
of the glass-forming liquids (empty magenta squares) is taken 
from [17fl . The thick black line is the gamma distribution with 
k = 30 and 6 — 0.175. The black dashed line is a Gaussian 
fit. The thin red line is the distribution of Voronoi structure 
of random points. The inset is the log-linear plot of the same 
data. 



Disorder and geometric constraint are two salient fea- 
tures of the jamming transition and the basic organizing 
themes underlying the jamming phase diagram. Both of 
them are essential in determining the shape of the uni- 
versal distribution, which can be seen as following. First, 
it is evident that for an ordered structure the area dis- 
tribution of Voronoi cells has zero variance and therefore 
should be a delta function (or delta functions). Second, 
we compare the distribution of our jamming system with 
that of totally random points without any geometric con- 
straint due to particle size (Fig.QJi). For random points, 
the peak of the distribution shifts to the left and the 
distribution has longer tail toward the right. 
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It is known that the Voronoi structure of glass-forming 
liquids also follows a universal distribution at different 
temperatures (T) ranging from liquid states to super- 
cooled states [17] . Even more striking, the universal dis- 
tribution of the Voronoi structure of glass-forming liq- 
uids is quantitatively similar to the universal distribu- 
tion of our 2D jamming granular system (Fig. [TJd) . It is 
likely that near the glass transition the interaction be- 
tween molecules of a supercooled liquid is dominated by 
its strong repulsive potential, therefore, like the ather- 
mal granular system, the geometric constraint of the 
molecules determines the local configurations and its dis- 
tribution. In addition, it is noted that deep inside the 
jammed phase the area distribution is invariant, which is 
consistent with previous study of the Voronoi structure 
of static granular packs at several discrete <j> [18j . Here, 
we show the invariant distribution with <j> varying con- 
tinuously and uniformly in a larger range. The universal 
distributions of Voronoi structure in glass- forming liquids 
at different T and in a granular system going through the 
jamming transition at different <f> corroborate the concept 
of the jamming phase diagram Qj. 

We also investigate the shape of Voronoi cells. In 2D, 
the geometric Euler constraint implies that average num- 
ber of edges of a cell should be 6 [H, ■ Cells with the 
number of edges other than 6 are usually referred to as 
"defects" 

(23423- 

We count the fractions of pentagons, 
hexagons and heptagons in the system (/s,/6 and f?) as it 
goes through the jamming transition. As shown in Fig. [2] 
fe increases while f$ and fj decrease: with increasing <f>, 
a pentagon-heptagon pair can anneal to two hexagons. It 
is interesting that the annealing process continues even 
at <fi clearly above the jamming point. The shape of the 
Voronoi cells is highly mobile due to tiny deformations 
of particles inside the jammed phase. Nevertheless, the 
annealing slows down deep inside the jammed phase and 
eventually stops. Due to the disordered nature of the 
jammed phase, a finite amount of the pentagon-heptagon 
pairs persist [20|. 

It has been proposed that the glass transition of 
binary-mixture molecular liquids can be characterized by 
the disappearance of certain kinds of "liquid-like" defects 
(heptagons enclosing small molecules and penta gon s en- 
closing large molecules) in Voronoi tessellation [2l|, ]22 1 . 
It is possible that in our granular system the "liquid-like" 
defects diminish during annealing and the remaining de- 
fects inside jammed phase are associated with "glass-like" 
defects. Due to the continuous distribution of the size of 
our particles and the difficulty to measure the accurate 
size of particles under water, we cannot directly verify 
the above scenario. However, a consistency check can 
be made by analyzing the average area of different types 
of cells. If heptagons enclosing small molecules, which 
should have smaller Voronoi cells in a condensed phase 
on average, diminish, then the average area of heptagons 
normalized by the total average area, should increase. 
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FIG. 2: Main plot: Number fraction of different type of 
cells with increasing <f>. Black squares for hexagons (/6), red 
circles for pentagons (/s) and blue triangles for heptagons 
(fr)- Left inset: Shape of cells in Voronoi tessellation of a 
pack at (j> = 0.86. Yellow cells are hexagons; blue cells are 
pentagons and red cells are heptagons. There is also one 
octagon (green) and one square (magenta). Note that most 
of pentagons and heptagons join into neighboring pairs and 
can anneal to hexagons later. Right inset: the average area 
of hexagons, pentagons and heptagons over the total average 
area of cells. The vertical dashed lines in both plots indicate 
the jamming point obtained from the structural signature of 
the jamming transition explained in the text. 



Similarly, if pentagons enclosing large molecules disap- 
pear, then the average area of pentagons over the total 
average area should decrease. We measured the normal- 
ized average area of pentagons, hexagons and heptagons, 
(s 5 )/(s), (se)/(s) and (s 7 )/(s), as a function of 4> in our 
experiments (Fig. fright inset). First, we can see pen- 
tagons have smaller area on average ((ss) / (s) < 1) and 
heptagons have larger area on average ({sr}/(s) > I). 
The trends of (s5)/(s), (se)/(s) and (sr)/(s) are more 
interesting. Before jamming the trends depend on the ini- 
tial packing configuration of individual experiment, sim- 
ilar to the behavior of P(s) before jamming. However, 
close to and after the jamming transition their features 
are robust for all experiments: (sq) / (s) plateaus, but 
(s5)/(s) clearly decreases and (s 7 ) / (s) increases. This 
is consistent with the behavior of glass-forming liquids 
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Before we can explore the relation between the glass 
transition and the jamming transition further, we need 
to ask a simple but important question: where does an 
athermal system jam? Though the jamming point of a 
pack can be easily identified by its mechanical properties 
|3j, experimentally it is usually hard to acquire further 
information of a pack except the position of all its parti- 
cles. Can we detect the jamming transition of a granular 
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FIG. 3: Pair correlations of a pack at different <j>. r is nor- 
malized by the original size of pearls do- The height of the 
first peak, g\, is indicated. Inset: The height of the first peak 
gi(4>) (blue squares) and the force measured at the bound- 
ary of the cell F{<f>) (red circles). The vertical dashed line 
indicates the peak of gi(<f)), <j>* ■ The horizontal dashed line 
indicates the zero force. 



pack from its structure? Despite the interesting features 
shown above, the Voronoi tessellation does not provide a 
unambiguous standard of the jamming transition. Both 
the area distribution and the annealing of defects show 
only quantitative changes before and after the transi- 
tion. We shall show next that, for a pack of granular 
particles with purely repulsive interactions, its pair cor- 
relation provides a structural signature of the jamming 
transition. 

We plot the pair correlation functions, g(r), of a pack 
at different <f> through the jamming transition (Fig. [3]). At 
each <j), g(r) has an oscillating shape, typical for amor- 
phous packs of spherical particles. Peaks in g(r) indicate 
the layering structure in the pack, the positions of which 
shift to larger r as <fr increases due to the enlargement of 
particles. A structural signature of the jamming transi- 
tion can be obtained by measuring <?i, the height of the 
first peak of g(r). gi((f>) shows a non-monotonic behavior 
with a peak at 0* = 0.84 ± 0.02 (Fig. [3] inset). We also 
measured the force applied by the pack on the boundary 
of the cell (Fig. [3] inset). The jamming point of the sys- 
tem is marked by the emergence of a non-zero force on 
the boundary, which is coincident with the peak of g\ (4>) ■ 
Hence, gi(<fi) provides a clear structural signature of the 
jamming transition and <fi* is the jamming point. 

The non-monotonic behavior of g\{4>) can be seen as 
the vestige of the structural singularity uncovered in the 
, amming transition of monodisperse frictionless particles 
ij. As a system evolves from unjammed state toward 



<j>*, particles are pushed closer to each other. The distri- 
bution of the distance between neighboring particles be- 
comes narrower and g\ therefore increases. For monodis- 
perse frictionless particles, due to the finite range purely 
repulsive interaction, at the jamming point the distances 
between contact particles are all exactly one particle di- 
ameter. Hence, gi, representing the probability to find a 
neighboring particle at a fixed distance, diverges. How- 
ever, it should be noted that for a system consisting of 
particles with attractive interactions the onset of com- 
plete connectivity of all particles and the rigidity of the 
pack need not occur at the same packing fraction. Above 
4>* , due to the disordered nature of the pack, particles will 
experience different compression and deform differently 
depending on their local environment. The distances be- 
tween contact particles spread out again. Quantitatively, 
<7i -w ~ z above jamming, where w is the left-hand width 
of the first peak of g(r) and z is the average coordination 
number of particles [4] . Since z increases slowly, it can be 
treated approximately as a constant near (f>* . Hence, as w 
increases from zero above </>*, g\ decreases. Apparently, 
the polydispersity of our particles depresses the diver- 
gence of g\ . Even at (p* , the distribution of the distances 
between contact particles still has a finite width. Never- 
theless, a non-monotonic peak is preserved as a vestige 
of the singularity. A simple estimation can be made on 
g\{<j>*) with polydisperse particles. Due to the difficulty 
to accurately measure the size of swelled soft particles, 
we estimate the polydispersity of swelled particles at <f>* 
to be around 5 to 8 percent, which is slightly smaller than 
that of dry particles. If we assume that there is no corre- 
lation between the radius of one particle and the radius 
of its neighbors, w should be roughly equal to the poly- 
dispersity of particles, i.e. w — 5% ~ 8% at <j)* . Using 
the relation gi(w) obtained from the simulations (Fig. 2 
of Ref.JH), we found that 2.8 < g\{4>*) < 4.5, which is 
consistent with our measurement (Figj3] inset). A similar 
thermal vestige of the structural singularity has also been 
observed in a colloidal system of finite T 14(, where the 
structural singularity is depressed by the thermal motion 
instead. 

By analyzing the Voronoi tessellation and the pair cor- 
relation functions of a granular pack with increasing 0, 
we have investigated the packing structures of the jam- 
ming transition. We have shown that the local configu- 
rations of a jamming granular pack are similar to those 
of glass-forming liquids, thereby demonstrating the anal- 
ogy between the jamming transition and the glass transi- 
tion in the spirit of the jamming phase diagram. Further 
study of the local structure of other jamming systems, 
such as glassy colloidal suspensions or jammed foams, 
may shed more light on the universal properties of the 
jamming structure illustrated here. Moreover, we have 
shown that a structural signature of the jamming tran- 
sition can be obtained from the first peak of pair cor- 
relation functions. Our study shows that, despite the 
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presence of realistic conditions such as frictional contacts 
and non-perfect spherical particle shapes, the structural 
signature observed in the ideal jamming system still per- 
sists in actual experiments. 
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